Discovery of a Strongly Lensed Galaxy at z = 3.9 behind a z = 0.83 
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Abstract 

We report the discovery and spectroscopic confirmation of three gravitationally-lensed images of a 
galaxy at z = 3.9 in the background of a distant, rich cluster of galaxies at z = 0.83, on the basis of 
observations with Faint Object Camera And Spectrograph (FOCAS) on the Subaru telescope. We construct 
a simple lens model of the cluster mass distribution based on Jee et al.'s weak lensing mass estimates from 
deep, high-resolution images by Advanced Camera for Surveys (ACS) on the Hubble Space Telescope. The 
lens model can account simultaneously for the observed image configuration and the flux ratio of the closer 
pair located close to the critical curve. The parities of the three images are also consistent with the lensing 
hypothesis. Since this galaxy is apparently bright (i^B ~ 23.7) for its redshift due to the magnification, 
it serves as a good high redshift target on which we can make extensive and detailed studies based on 
multi- wavelength imaging and spectroscopy. 
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1. Introduction 

Strong gravitational lensing due to the dense cores of 
clusters of galaxies leads to the formation of multiple im- 
ages and giant luminous arcs of background galaxies. The 
lensing observables such as image positions and flux-ratios 
in turn provide strong constraints on the underlying grav- 
itational potential of the lensing clusters. The strongly 
lensed images of background galaxies have therefore been 
powerful observational tools to probe the mass distribu- 
tion in cluster cores (Mcllier et al. 1993; Kncib et al. 1993, 
1996, 2003; Hattori, Makino, & Kneib 1999; Molikawa & 
Hattori 1999; Broadhurst et al. 2005). 

Exploring high- z galaxies is extremely important for un- 
derstanding unsolved physical properties of forming galax- 
ies as well as understanding the cosmic reionization his- 
tory. The bottleneck is that such high-z forming galaxies 
are too faint to be observed with present-day telescopes 
due to the (1 + z) 4 -dependence of the surface-brightness 
dimming effect. The most promising method is to utilize 
the gravitational magnification by clusters as a biggest 
natural telescope - Gravitational Lens Telescope (Ellis et 
al. 2001; Hu et al. 2002; Kneib et al. 2004). Furthermore, 
detailed study of high redshift galaxies can only be made 
for strongly lensed hence magnified galaxies (e.g., Pettini 

* Based in part on data collected at the Subaru Telescope, which 
is operated by the National Astronomical Observatory of Japan. 



et al. 2002 on cB58). 

We have discovered a high-z lensed galaxy during the 
course of spectroscopic follow-up of the RXJ0152.7— 1357 
cluster at z = 0.83 as a part of our distant cluster project 
on the Subaru called PISCES (Kodama et al. 2005). 
In this paper, we present the spectroscopic and photo- 
metric properties of the newly discovered triple lensed 
images based on the observational data taken with the 
ground-based Subaru telescope as well as on the data 
with the high-resolution archival Advanced Camera for 
Surveys (ACS) on the Hubble Space Telescope (HST). 
A gravitational lens model of the cluster mass distribu- 
tion is developed for the three lensed images. Throughout 
the paper we assume a A-dominated flat cosmology with 
TO = 0.27, Oa = 0.73, and h = 0.7, and we use magnitudes 
given in the AB system. 

2. Identification of the Lensed Objects 

The lensed objects were first found as those having very 
similar colors in VRi'z' based on the optical imaging with 
Suprime-Cam on the Subaru telescope (see Kodama et al. 
2005 for details of the optical data). The photometric 
properties of the objects are summarized in Table 1. All 
the three objects show strong continuum breaks in the 
V — R color (~ 0.9) with a relatively flat SED at longer 
wavelength regime. These photometric properties suggest 
that their redshifts are high (z ~ 4) in analogy to Lyman 
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Fig. 1. The optical spectra of three multiple images 
(7i, I2, I3) obtained with FOCAS on the Subaru telescope. 
The zero-flux levels are indicated with horizontal dotted lines. 

break galaxies. 

We have conducted spectroscopic follow-up observa- 
tions during 11-14 October 2004 with FOCAS (Kashikawa 
ct al. 2002) in the MOS mode. We used a 300 lines mm -1 
grating blazed at 5500 A with the order-cut filter SY47. 
The wavelength coverage was between 4700A and 9400 A 
with a pixel resolution of 1.40 A pixel -1 . A slit width was 
set to 0"8, which gave a resolution of A/AA ~ 500. We 
obtained four 1800s exposures for each object nodded by 
~ 0."6 for each shot. The object I 2 was slipped out of the 
slit for one shot due to the nodding, and the net exposure 
on this object is 3 x 1800s. Data reduction was performed 
in a standard manner using IRAF. The objects I\ and 1 2 
were very close to the spatial edge of the slits, and hence 
we will not discuss relative fluxes of these objects from the 
spectra. Figure 1 shows the obtained spectra of the the 
three images under concern. All the three show a strong 
emission line around A b s = 5992A. The emission line pro- 
file for each object was fitted with a Gaussian. The cen- 
tral wavelengths of the lines are 5992.41A, 5991.95A, and 
5992. llA for I\, I 2 , and I3, respectively, with a statistical 
error of ~ 0.3A. Therefore, the objects lie at the same 
redshift within the errors. The objects I\ and -Z3 show 
a clear continuum break shortward of this line, although 
no such break is seen in I 2 due to the low S/N. No other 
lines are observed in the spectra. These objects are ei- 
ther (1) a gravitationally lensed Lya emitter at z = 3.928 
or (2) foreground [Oil] emitters at z = 0.608. The lat- 
ter possibility is rejected, however. We observe no [OIII] 
lines at the expected wavelength as shown in Fig. 1. It 
is highly unlikely that galaxies show such a strong [Oil] 
emission without any detectable [OIII] emission. The pho- 
tometric properties of the objects are inconsistent with a 
z = 0.608 galaxy either. Their V — R and R — i' colors are 
too red for their i' — z 1 colors to be star forming galaxies 
at z ~ 0.6. We therefore conclude that these objects are 
Lya emitters at z = 3.928. Note that the rms error in the 
redshifts is <r(z) ~ 3 x 10~ 4 . In summary, these objects 
show very similar photometric properties and they lie at 



Fig. 2. Zoom-in view of a 20" X 20" region around the three 
gravitationally-lcnsed objects (ii, ^2,^3) in the HST/ACS 
composite color image constructed using F675W, F775W, and 
F850LP filters. The stick at the bottom right has a 5" length. 
The red open rectangles indicate the slit locations for the 
Subaru/FOCAS multi object spectroscopy of (Ii, I2, 13). The 
inset plots show the surface brightness contours of {Ii,l2,l3) 
from the F675W-band data. The coordinates are in arcsec- 
onds. The lowest contour and the contour interval are both 
at a 3<t rms noise level of 25.06 mag arcsec - 2 . 

exactly the same redshift - these facts strongly suggest 
that they are originally a single galaxy strongly lensed by 
the foreground cluster at z — 0.83. 

We also obtained the "on-the-fly" processed HST/ACS 
images of the field from the HST archive (HST GTO 
Proposal 9290: Ford, H.). The cluster was observed in 
a 2 x 2 mosaic pattern with an overlap of ~ 50" regions in 
F675W, F775W, and F850LP filters. The total exposure 
per pointing is about 4.8 ks for each filter. We show in 
Fig. 2 the HST image of the lensed objects and the Subaru 
image in Fig. 3. The objects all apparently look like the 
edge-on disk, which is probably due to the effect of strong 
lensing. 

3. Strong Gravitational Lensing 

3.1. Lens Model 

The RXJ0152 gravitational lens system shows the col- 
incar, three multiple images I 2) I^) of a source galaxy 
at z — 3.9, from which one expects that the lensed images 
could be associated with a cusp caustic produced by an 
elliptical potential (e.g., Schneider, Ehlers, & Falco 1992). 
In constructing a lens model of the cluster mass distribu- 
tion, we closely follow weak lensing mass estimates by Jee 
et al. (2005) based on the HST/ACS observations with 
superb angular-resolution and photometric-depth. We 
adopt Jee et al.'s parameterized singular isothermal ellip- 
soid model (SIE: Kormann et al. 1994; King & Schneider 
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Table 1. Photometric properties of three gravitationally lensed images from Subaru/Suprime-Cam and HST/ACS data* 



Object 


R.A. 


Dec. 


»t 


V-R 


R-i' 


i'-z' 






(J2000) 


(J2000) 


(AB-mag) 


(AB-mag) 


(AB-mag) 


(AB-mag) 


(") 


h 


01 52 45.32 


-13 57 08.6 


23.70 ±0.04 


+0.85 ±0.01 


-0.03 ±0.01 


-0.07 ±0.03 


0.141 


h 


01 52 45.14 


-13 57 04.9 


24.05 ±0.04 


±0.90 ±0.02 


±0.08 ±0.02 


+0.18 ±0.03 


0.125 


h 


01 52 44.41 


-13 56 56.7 


23.89 ±0.04 


±0.88 ±0.01 


-0.06 ±0.02 


+0.02 + 0.03 


0.127 



* The total magnitudes 'MAG_AUTO' (SExtractor) are used for 2^b- Photometric apertures of l" are used for color measurements. 
% Half-light radii (FLUX .RADIUS by SExtractor) measured in the HST/ACS F675W-band image 
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Fig. 3. VRi' color image of the central 90" X 90" region 
of RXJ0152 (z = 0.83) obtained with Subaru/Suprime-Cam. 
North is top and East is left. Also overlayed are the caustic 
(thick solid lines) and the critical curves (thin dashed lines) 
predicted by our best fitting model. The cross represents 
the predicted position of the source galaxy at z = 3.9. The 
open squares represent the predicted positions of the multiple 
images (h,h,h)- 

2001) centered between the two central brightest cluster 
galaxies, for describing the elongated global mass distribu- 
tion. The SIE requires three model parameters: Einstein 
radius #e, minor-to-major axis ratio /(= b/a), and orien- 
tation angle a. Jee et al. obtained the best-fit param- 
eters of E = 7"19 ± 0^72 (z s = 1.3) and / = 0.36 ± 0.1. 
This Einstein angle corresponds to the 1-D velocity dis- 
persion of tr,^ IE = 940 ± 168kms _1 . For the source redshift 
of z s = 3.9, this angle is scaled to 9e = 15" 4. The cluster 
contains several mass clumps associated with the cluster 
galaxy concentrations. As the three images are located in 
between the cluster center and the North-East substruc- 
ture (mass clump B of Jee et al.), we take this clump into 
account in our lens model in the following way: We as- 
sume the clump mass of M = 3.4 x 1O 13 M estimated by 
Jee et al. is the virial mass of the clump. We then con- 
vert the mass into the 1-D velocity dispersion of an SIS 
halo (of IS oc M 1 / 3 : see Bartelmann, King, & Schneider 



2001). Accordingly, our model consists of one SIE as the 
global cluster halo and one SIS as the North-East sub- 
structure. We treat the SIE parameters (#e, /, a) and 
the (un-lensed) source position 0s as free parameters of 
our lens model. 

We then constrain the 5 model parameters via x 2 -fitting 
to the observation data: the positions and the fluxes of the 
three images. In the model-fitting, wc use as constraints 
the thee image positions Oi k (k = 1,2,3) relative to the 
center of the SIE halo; we further include the flux ratio 
I2 1 1\ of the closer pair as it gives a strong constraint on 
the location of the critical curve. The x 2 is then given as 



E 

fc=i 



, [{hih)-{h/h) M \ 



(i) 



where quantities with superscript M denote the model pre- 
dictions, cr pos is the positional uncertainty of the multi- 
ple images, er pos wfll'l, and uj 2 /j 1 is the rms error of the 
flux ratio I2 /h ■ The image positions 61 of a source with 
the position 83 are obtained by solving the lens equation 
8 s = 0i- Vip(8i) (Schneider et al. 1992) with the effec- 
tive lensing potential ip = tpsm + V'sis m sum of an SIE 
and an SIS potential. 

3.2. Model Predictions v.s. Observations 

We show in Fig. 3 the derived best-fit model of the 
RXJ0152 lens system overlaied on the Subaru image. 
The minimum x 2 is Xmin/dof = 2.03/2 for 2 degrees 
of freedom (dof). We quote the best-fit SIE param- 

14 // 3 +0.4( + l.l 



eters are Of 
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a = 47.°7±0.2(±0.4) measured East of North, with the 
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la (2a) errors estimated from A^ 2 = X 2 ~ X 
the 5 parameter space. 

The predicted source position is inside and close to a 
cusp caustic produced by the main cluster potential as 
expected. For such an image configuration, the image I2 
is predicted to have opposite parity as I\ and ^3 . In partic- 
ular, the closest images (/i,/2) are mirror images of each 
other with respect to the critical curve (see Fig. 3). The 
high-resolution ACS images reveal more detailed struc- 
tures of the lensed images that were not spatially resolved 
in the Subaru images. Figure 2 clearly shows that the ob- 
served image parities are consistent with the lensing hy- 
pothesis: each image contains a bright core and extended 
emission with lower surface brightness, and a clear mirror 
symmetry between I\ and I2 is seen in the HST/ACS data. 
The predicted amplifications using the best-fit model are 
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/Xi = 7.2, fa = 5.8, and /i3 = 3.0 for ii,/2, and ^3, respec- 
tively. The total amplification factor over the three images 
is hence w 16. Accordingly, the predicted magnification 
ratios are J2/J1 = 0.80 and I3/I1 = 0.42, whereas the flux 
ratios from the Subaru i'-band data are I2/I1 = 0.82 ±0.04 
and I3/I1 = 0.89 ±0.05. Our simple lens model thus re- 
produces correctly the observed flux ratio of I2/I1 while 
it fails for I3/I1; more than a factor 2 smaller than the 
observed flux ratio. In this study, we will not assess the 
significance of the difference between the observed and 
predicted flux-ratios for ^3 / 1\ as it requires more detailed 
lens modeling for constraining all possible substructures 
(cf. Kncib ct al. 1996). It is expected that the image 
73 could be easily magnified by local perturbers since its 
angular size (see Table 1) is sufficiently small compared 
with that of the lensing potential of possible perturbers 
(Chiba 2002). 

4. Intrinsic Properties of the Source Galaxy 

The derived lens model allows us to inspect the intrin- 
sic properties of the un-lensed source from the observed 
images. In what follows, we only refer to the images 
I\ and I2 for which our lens model can accurately re- 
produce the observations. From the predicted magnifica- 
tions, the magnitude of the un-lensed source is estimated 
as i' AB ~ 25.9. This corresponds to a UV-luminosity of 
F(1560A) = 6.3 x 10 40 erg/s/A or M(1560A) = -20.2 mag. 
From this, we estimate a star formation rate (SFR) of 
the source galaxy to be 6MQyr _1 without dust extinction 
correction (Madau et al. 1998). The observed Lya flux 
indicates SFR~ 3M Q yr _1 (Kennicutt 1998; Brocklehurst 
1971), which is lower than SFR(1560A). This is suggestive 
of dust extinction. Further, the superb HST/ACS angular 
resolution together with our lens model gives us size infor- 
mation of the source galaxy. The half-light radii of the 
lensed images are listed in Table 1 (FLUX_RADIUS given 
by SExtractor: Bertin & Arnouts 1996). Note that even 
for the ACS images the FWHM of PSF is comparable to 
that of the lensed images. Taking this into account, we 
give a constraint on the intrinsic source size as < 0"02 
(>3a <0."08; see Fig. 2), or < 150 pc (<600 pc) in physical 
scales, assuming a simple one-dimensional amplification 
near the critical curve (see Ellis et al. 2001 for a similar 
case). 

5. Summary 

We report the discovery and spectroscopic confirmation 
of three gravitationally lensed images of a distant galaxy 
in the background of the rich cluster RXJ0152.7— 1357 
(z = 0.83). The similarity of the photometric and spec- 
troscopic properties of the three images from Subaru ob- 
servations provides strong evidence that they are gravita- 
tionally lensed images of a Lya emitter at z = 3.9. A lens 
model was constructed on the basis of Jee et al.'s weak 
lensing mass estimates from the HST/ACS data. Our 
simple lens model can account simultaneously for the ob- 
served image configuration and the flux ratio of the closer 



pair located close to the critical curve. The image pari- 
ties in the ACS image are also consistent with the lensing 
hypothesis. 

To further investigate the properties of this source 
galaxy at z = 3.9, we need longer wavelength data, such 
as A'-band which samples its rest-frame B-band light or 
Spitzer/ Astro- F bands at 3-8/im which sample the light 
from underground old populations related to the stellar 
mass of this galaxy. A more detailed modeling of the 
lens mass distribution is important, as a well-constrained 
cluster mass model in conjunction with different sets of 
multiple images with suitable redshifts serves as a pow- 
erful probe of the geometry of the universe (Futamase & 
Yoshida 2001; Golse, Kneib, & Soucail 2002). 
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